Anomaly in Nonlinear Magnetoelectric Response of YbMnOs 
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We observe a seemingly complex magnetic field dependence of dielectric constant of hexagonal 
YbMnOa near the spin ordering temperature. After rescaling, the data taken at different temper- 
atures and magnetic fields collapse on a single curve describing the sharp anomaly in nonlinear 
magnetoelectric response at the magnetic transition. We show that this anomaly is a result of the 
competition between two magnetic phases. The scaling and the shape of the anomaly are explained 
using the phenomenological Landau description of the competing phases in hexagonal manganites. 
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The recent interest in multiferroic materials was trig- 
gered by the discovery of the giant magnetocapacitance 
(MC) and magnetically-induced rotations of electric po- 
larization in orthorhombic rare earth manganites 
This multiferroic behavior is rooted in magnetic frustra- 
tion, which gives rise to non-centrosymmetric spin or- 
derings that induce electric polarization [4|. Further- 
more, the presence of competing spin states in these 
frustrated magnets results in a strong sensitivity of 
the magnetically-induced electric polarization to applied 
magnetic fields. In this respect multiferroics are sim- 
ilar to colossal magetoresistance manganites and high- 
temperature superconductors Q. 

In this Letter we study effects of critical magnetic fluc- 
tuations and the competition between different magnetic 
states on the non-linear magnetoelectric response of the 
hexagonal YbMnOa by measuring the magnetic field and 
temperature dependence of its dielectric constant. Ferro- 
electricity in hexagonal manganites i?Mn03 (i? = Ho-Lu, 
Y) appears well above the magnetic transition and is of 
nonmagnetic origin: An electric dipole moment along the 
c axis is spontaneously induced by tilts of manganese- 
oxygen bipyramids and buckling of rare earth-oxygen 
planes at Tc > 600K 0,0, H,!!, while the ordering of Mn 
spins occurs at a much lower temperature T/v < 120K. 
However, the spin ordering in hexagonal manganites re- 
sults in a surprisingly strong lattice relaxation, which 
affects the spontaneous electric polarization [lo| . 

The Mn ions in hexagonal manganites form well- 
separated triangular layers parallel to the ab plane with 
antiferromagnetic exchange interactions between nearest- 
neighbor spins 11], which makes the Mn spin subsystem 
low-dimensional and frustrated and results in enhanced 
spin fluctuations observed well above T/v [l2|- Frustra- 



tion and rare earth magnetism are responsible for a rich 
variety of magnetic phases observed at low temperatures 
and in applied magnetic fields [l3|- Due to magneto- 
electric coupling each magnetic transition gives rise to a 
singularity of the dielectric constant [1, Q, [3, , which 
is more pronounced than the corresponding singularity in 



magnetic susceptibility. 

We find that close to the Neel temperature Tat— 81K 
the MC of YbMnOs measured as a function of magnetic 
field and temperature obeys a scaling behavior and has 
a very sharp anomaly. The detailed comparison with re- 
sults of model calculations led us to a conclusion that the 
effect of magnetic fluctuations is completely overshad- 
owed by the magnetic field dependence originating from 
the competition between two antiferromagnetic states, 
one of which is weakly ferromagnetic. Using a mean field 
Landau expansion of free energy in powers of two com- 
peting order parameters, we reproduce the shape of the 
anomaly as well as the changes in the behavior of MC 
observed in the wide range of magnetic fields and tem- 
peratures. 

Polycrystalline samples of YbMnOs were prepared by 
solid state synthesis. A single crystal was grown from 
this powder by the floating zone technique. Magneti- 
zation M(T) of the samples was measured by a Squid 
magnetometer (MPMS7 Quantum Design) using a field 
of 0.5T. Field dependence of the magnetization was mea- 
sured up to 5T. Capacitance of the samples was measured 
in an commercial system (PPMS Quantum Design) us- 
ing a home-made insert and a Andeen-Hagerling 2500A 
capacitance bridge operating at a fixed measurement fre- 
quency of 1 kHz as well as using an Agilent 4284A LCR 
meter up to IMHz. Electrical contacts were made using 
Ag epoxy. 

The temperature dependence of the capacitance C (T) 
proportional to the in-plane dielectric constant Sa is 
shown as an inset of Fig. [T^. Below the Neel tem- 
perature Tj\f — 82K, the capacitance is somewhat sup- 
pressed by the emergence of magnetic order [l^. The 
MC, where H is magnetic field along the c 

axis, for a set of temperatures between 76. 5K and 95K 
is shown in Fig. [1^. In this small temperature interval 
around T/v the behavior changes dramatically: at 76. 5K 
only a positive curvature is observed. With increasing 
temperature a high-field downturn appears, and at 80K 
only a negative curvature can be observed, which changes 
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back to positive above 90K. 

This unusual behavior is a consequence of the fact that 
by varying temperature and magnetic field we force the 
system to pass through a magnetic transition. The criti- 
cal behavior becomes apparent when we plot '^g^^o)//!"'' 
vs temperature [see Fig. mb)]. The procedure to evalu- 
ate AC at fixed magnetic fields versus temperature i.e. 
replotting the rescaled changes of dielectric constant in 
magnetic field versus temperature, effectively reveals the 
magnetic field dependence of C. The strong temperature 
dependence of C masks the magnetic field dependence 
when C(T) is measured at fixed magnetic fields. The 
data taken at various T and H remarkably fall onto a 
single curve with a very sharp anomaly at T^r where the 
temperature derivative of MC becomes large and posi- 
tive while its magnitude shows an almost discontinuous 
jump from a positive to a negative value. The observed 
scaling behavior of MC can be understood in terms of 
the anomalous nonlinear magnetoelectric response, de- 
scribed by the term n{T) (^E^ + i?^) in free energy, 
where k{T) has a singularity at the magnetic transition 
temperature. The scaling implies that the dependence of 
magnetic susceptibility for H\\c on electric field E\\a has 
the same anomaly at T^. 

The shape of the MC anomaly in YbMnOs is un- 
usual. The nonlinear magnetoelectric response in an- 
tiferromagnets usually originates from the ubiquitous 
fourth-order coupling of the electric polarization P to 
the magnetic order parameter L, /,nc = ^P^L'^- It re- 
sults in a correction to the bare dielectric susceptibility 



Tjv 



< 0, which 



Xo, Sx = -gxl {L) oc (-r)2P, for r 
accounts for the observed dielectric constant anomaly be- 
low T/v [see the inset of Fig. [Ifa)] . The magnetic field de- 
pendence of Neel temperature, Tjv {H) « Tjv (0) — ATJ^, 
gives rise to a discontinuity of MC at Tn and its anoma- 
lous behavior below Tn, is roughly consistent with our 
data. However, the most prominent feature of the ob- 
served anomaly - the long negative tail for T > T/v [see 
Fig. [TJb)] - cannot be explained in this way. 

This tail may result from magnetic fluctuations that 
become critical close to Neel temperature. The two 
lowest-order self-energy diagrams describing contribu- 
tions of magnetic fluctuations to dielectric susceptibil- 
ity are shown in Fig. [21 The lowest-order term is given 

by i^x''^'' = —5X0 {{^^) ^ W'^^ °^ T^^'^j where a is 
the exponent describing the critical behavior of mag- 
netic specific heat [T^- The corresponding singular- 
ity in MC oc 5Asign(T) Irp" (g, A > for YbMnOg), 
is positive for T > Tjv, in disagreement with our data. 
The second-order fluctuational correction [see Fig. [^b], 
^^(2) ^ ^L9M^Jd3x\{L^^)L^0))-{L^f] , results 

in the MC anomaly oc — Asign(r) |r| ^ which has the 
shape similar to that observed experimentally. However, 
this term is only nonzero in a ferroelectric state with a 




FIG. 1: (Color online) (a) Magnetic field dependence of MC 
of YbMnOa single crystal at constant temperatures near Tat. 
The temperature dependence of capacitance is added as an 
inset, (b) Temperature dependence of MC at constant mag- 
netic fields. The inset shows the shift of Tjv in magnetic field 
obtained from magnetic susceptibility measurements. Electric 
field is parallel to the ab plane, while H\\c. 



(a) 





FIG. 2: (Color online) The self-energy diagrams of the first 
(a) and second (b) order in the coupling constant g, describing 
contributions of magnetic fluctuations (wavy lines) to dielec- 
tric susceptibility. 



spontaneous polarization Pq ^ and, therefore, should 
be visible in Sc, while we observe an anomaly only in ea,b- 

From this we conclude that the observed anomaly is 
unrelated to magnetic fluctuations and originates from 
a different physics. Below we show that the shape and 
scaling behavior of MC can be explained within a mean 
field theory by the competition between antiferromag- 
netic and weakly ferromagnetic state. 

Hexagonal manganites show a number of magnetic 
phases with the 120°-angle between Mn spins in triangu- 
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FIG. 3: (Color online) The magnetic phase diagram of the 
model Eq.ll} for T^°'' > The hlac region is the B2 phase 

with some admixture of the A2 phase. The critical point 
(CP) separates the first-order transition (solid) line from the 
second-order transition (dashed) line. 



lar ab layers [l^ . These phases differ by orientation of the 
spins with respect to the crystallographic axes and spins 
in neighboring' layers as well as by the ordering of rare 
earth spins [19, JQi, [21], H^l- The magnetic phase diagram 
of YbMnOa studied by a variety of different experimen- 
tal techniques includes the low-field B2 phase (magnetic 
space group PGscm) and the high-field A2-phase (mag- 
netic space group PBaCin) [H, H^l . The symmetry of the 
latter state allows for a net magnetization in the c di- 
rection (largely due to the rare earth spins and therefore 
small near the Mn spin-ordering temperature), which is 
why the A2 phase is stabilized by H\\c. 

The competition between the A2 and B2 phases was 
discussed in Ref. 24] using the phenomenological free 
energy expansion in two order parameters: 



/ 



E 

7=A,B 
d 



T 



rp(0) 

7 



r2 t2 



H (j)Lp 



3 



"77-4 



(1) 



where La{Lb) is the order parameter describing the 
A2(B2) phase and H is magnetic field along the c axis. 
The linear coupling of Xa to H corresponds to the spon- 
taneous magnetization present in the A2 phase. 



The typical phase diagram for r]^°' > T|"^ (when the 
B2 phase is energetically more favorable than the A2 
phase at zero field) is shown in Fig. [H Due to the linear 
coupling between H and La, the latter order parameter 
is nonzero for an arbitrarily weak magnetic field, so that 
ioi H ^ {) the transition occurs between the A2 phase 
and the B2 phase with some admixture of the A2 phase. 

The MC for this model, shown in Fig. Hljc), is calcu- 
lated by adding to the free energy Eq.([T|) the terms de- 
scribing the coupling of the magnetic order parameters 



.(0) 





(a) 




. "-A 






T 





La 




(b) 


H 



I 



I 

X 





(c) 






T 




^ 

(d) 



H 



FIG. 4: (Color online) The temperature dependence of order 
parameters La (dashed blue line) and Lb (solid red line) near 
the second-order (a) and first-order (b) transitions. Panel 
(c) shows the anomaly in rescaled MC at the second-order 
transition temeprature; the two curves calculated for different 
values of magnetic field are almost indistinguishable. The MC 
anomaly at the first-order transition is shown in panel (d). 



to the in-plane electric polarization and the dielectric re- 
sponse of the nonmagnetic state, 
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PE, (2) 



has the same shape as the one observed in YbMnOa and, 
for weak fields, obeys the observed scaling. This behav- 
ior can be understood by noting that the main contri- 
bution to the magnetic field dependence of the dielec- 
tric susceptibility conies from La, which is linearly cou- 
pled to H . This field-induced order parameter grows as 
T approaches the T/v from above [see Fig. HJ^a)], which 
gives rise to the 'high-temperature' negative MC tail, 
as oc —L\. In the weak- field regime La oc so 
that ^'^'^'^^^^^''-'^''^^ is approximately field-independent, 
which explains the observed scaling [2a |. 

As the magnetic field increases, the character of the 
transition in the two-parameter model changes: in low 
fields the transition is of second order (red dashed line 
in Fig. [3]), while in high fields and low temperatures 
it becomes a first-order transition (blue line) [see also 
Figs. IDJa) and (b)]. The first- and second-order transi- 
tion lines are separated by the critical point. This change 
in the nature of the transition is also clearly seen in the 
experiments by comparing the field-dependence of MC 
at low temperatures [see Fig. EJa)] to that at high tem- 
peratures [see Fig.[T{a)]. At 2K the MC shows a distinct 
cusp at the first-order transition, which is well repro- 
duced within our model [see FigHJd)]. In YbMnOa the 
changes in the order of the transition are made more dra- 
matic by the fact that at low temperatures and high mag- 
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FIG. 5: (Color online) (a) MC of YbMnOa at T < Tjv for 
i5||a6-plane and il\c\ (b) Field dependence of magnetization 
at 2 and 5K. 



netic fields magnetic response is dominated by Yb spins, 
which below 3.8K order ferrimagnetically [2^. This leads 
to a strong decrease of the critical magnetic field at low 
temperatures (ilr_ ~ 3T at 2K) and gives rise to the sharp 
discontinuity in magnetization [see Fig. [HJb)] [26| . 

Finally, since the shape of the MC anomaly may be af- 
fected by the motion of antiferromagnetic domain walls 
separating different antiferromagnetic phases, we mea- 
sured the frequency dependence of the dielectric con- 
stant. We found no frequency dependence in the interval 
from IkHz to IMHz, suggesting that the contribution of 
the domain walls in YbMnOa is negligible. 

To summarize, we observed a sharp anomaly in nonlin- 
ear magnetoelectric susceptibility of the hexagonal rare 
earth manganite YbMnOa at the Neel temperature. We 
discussed theoretically possible sources of the anomaly 
and showed that it results from the competition between 
two antiferromagnetically ordered states of YbMnOs , one 
of which has a small spontaneous magnetic moment. 
Even though this weakly ferromagnetic phase becomes 
the ground state only in rather high magnetic fields or 
at very low temperatures, its admixture to the non- 
ferromagnetic phase determines the shape of the MC 
anomaly along the whole critical line of magnetic phase 
transitions. The competition between different magnetic 
phases, some of which may have weak ferromagnetic mo- 
ment is very common for frustrated magnets: similar 
phase diagram were found for hexagonal HoMnOs, which 
shows four competing states [3| and for Ni3V208 [27t . 
Thus many other systems should show an anomaly in 
nonlinear magnetoelectric response, although its shape. 



which depends on parameters of the Landau free energy, 
may vary from material to material. 
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